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The Cosmological Principle in it's various versions states that: (i) the Galaxy 
does not occupy a particular position, (ii) the Universe is homogeneous and 
isotropic. This statement does not agree with the recent astronomical obser- 
vations in the range z lower than 0.05 which are in agreement with a cellular 
structure of the Universe. Here we present a local analysis of the inhomo- 
geneity of the Universe. When z is greater than 0.05 our analysis cannot 
be applied because the astronomical sample of galaxies here processed is not 
complete. The two tools of the Poisson Voronoi tessellation (PVT) and the 
luminosity function for galaxies allow of building a new version of the local 
Cosmological Principle, 
keywords 

Cosmology: miscellaneous Cosmology: observations Cosmology: theory 



1 Introduction 

All the various theories which belong to cosmology are based on the "Cosmo- 
logical Principle" which is now presented in two versions. The first version 
asserts that "There is nothing special about our location in the Universe," 



see Ryden (2003); Liddle & Loveday (2009). The second version is "Viewed 



on a sufficiently large scale, the properties of the Universe are the same for 
all observers," see Keel (2007). Other versions which allow the formation of 



structures and incorporate probability distributions can be found in Schwarz 



(2009) 



These two statements are strictly connected with the concept of isotropy 
and 3D homogeneity of the Universe. In the following we will show that 

• The Local Super Cluster (LSC), the nearest part of the universe, has 
a dimension around 2 which is characteristic of a bi-dimensional space 
and therefore the concept of 3D homogeneity is violated by the obser- 
vations. 

• The concept of isotropy in the spatial distribution of galaxies is not 
confirmed due to the 3D cellular structure of galaxies. 

In order to demonstrate the previous two statements this paper is orga- 
nized as follows. Section [2] reviews the adopted cosmology. Section [3] reports 
a spherical cut of the 3D Voronoi diagram as well an introduction to the thick 



Voronoi faces. Section <L2 reports a spherical cut of an astronomical catalog. 
Section [5] reports the most used luminosity function for galaxies, a modern 
version of the Malmquist bias, and the measure of the spatial dimensions of 
the LSC. Section M reports a revised version of the "Cosmological Principle." 



2 Elementary Cosmology 



Starting from |Hubble| ( |1929 ) the suggested correlation between expansion 
velocity and distance is 

V = H D = cz , (1) 

where the Hubble constant is Ho = 100/i km s _1 Mpc -1 , with h = 1 when h 
is not specified, D is the distance in Mpc, c is the velocity of light, and z is 
the redshift. The redshift of galaxies is explained by the Doppler effect or by 
alternative theories such as the presence of a hot plasma with low density, see 



Brynjolfsson (2009). In a Euclidean, non-relativistic, and homogeneous uni- 

L G) 

verse, the flux of radiation, /, expressed in ^% units, where Lq represents 



2 



the luminosity of the sun, is 



AttD 2 ' 

where D is the distance of the galaxy expressed in Mpc, and 



D 



H 



(2) 



(3) 



3 Voronoi Diagrams 



The points of a tessellation in 3D are of four types, depending on how many 
nearest neighbors in ES, the ensemble of seeds, they have. The name seeds 
derive from their role in generating cells. Basically we have two kinds of 
seeds , Poissonian and non Poissonian which generate the Poissonian Voronoi 
tessellation (PVT) and the non Poissonian Voronoi tessellation (NPVT). The 
Poissonian seeds are generated independently on the X, Y and Z axis in 3D 
through a subroutine which returns a pseudo-random real number taken from 
a uniform distribution between and 1 . This is the case most studied and 



for practical purposes, the subroutine RAN2 was used, see Press et al. (1992 ) 



The non Poissonian seeds can be generated in an infinite number of different 
ways : some examples of NPVT are reported in Zaninetti (2009). 



A point with exactly one nearest neighbor's is in the interior of a cell, a 
point with two nearest neighbors is on the face between two cells, a point 
with three nearest neighbors is on an edge shared by three cells, and a point 
with four neighbors is a vertex where three cells meet. The Abell clusters 
and therefore the galaxies were originally inserted on the Voronoi vertexes, 
which are the points which share the same distance from four seeds, see 



van 



de Weygaert & Icke (1989). The target to reproduce in the previous analysis 
was the cluster-cluster correlation function. In order to explain the cellular 
structure of the universe the galaxies were later on collocated on the faces 
of a NPVT in order to reproduce the galaxy-galaxy correlation function and 
the cluster-cluster correlation function, see Zaninetti (1991). More detailed 



simulations of the various astronomical catalogs as well of the galaxy-galaxy 
correlation function were recently obtained inserting the galaxies on the faces 
of a PVT, 



sec 



Zaninetti (2006 2010). 



According to this trend we will assume that the galaxies are situated on 
the faces of the PVT. Some of the properties of the PVT may be deduced 
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from approximate arguments introducing the averaged radius of a polyhe- 
dron , R and the averaged diameter D = 2R. This theoretical quantity 
has its observational counterpart in the averaged diameter of the voids be- 

;he effective radius of the voids 



tween galaxies, D obs . A careful analysis of t 
between galaxies in SDSS DR7, see Section 
The averaged volume is 



4.1 



derives D obs = ^ Mpc 



V=~<kR 3 , (4) 



and the averaged number of faces is n which according to Okabe et al. ( 1992 ) 
is n = 15.35 . 

The averaged surface area of a polyhedron, S, is 

S = AttR 2 , (5) 

and the approximate of a face, A, is 

A=*^ . (6) 
n 

The approximate side of a face, f, is 



On assuming that our galaxy is at the center of an irregular face of a PVT, 
the galaxies on the other faces are comprised within a distance Dp 

D f m 1 - . (8) 

The cross sectional area of a PVT can also be visualized through a spherical 
cut characterized by a constant value of the distance from the center of the 
box, in this case expressed in z units, see Figure [TJ this intersection is called 
V s (2, 3) where the index s stands for sphere. The galaxies are thought to be 
situated on the faces of a PVT network and Figure [2] reports an example of 
a spherical cut. 

An astronomical slice consists of a pie diagram with polar angle given by 
the right ascension and polar radius given by the redshift. An example of a 
slice of the Sloan Digital Sky Survey (SDSS) which contains all the galaxies 
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Spherical VP, elements = 5865 . at 2=0.0133 




Galactic Longitude, 1 



Figure 1: The Voronoi diagram V s (2,3) in the Hammer-Aitoff projection 
0.0133 < z < 0.0144 or 56.9 Mpc < D < 61.67 Mpc when #=18.16 Mpc. 



Spherical VP, elements^ 693 , at 2=0.0133 




Ga lactic Longitude, I 



Figure 2: The Voronoi diagram V s (2,3) in the Hammer-Aitoff projection 
when 693 galaxies are extracted from the network of Figure [TJ 
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within an opening angle in declination ,ip , of 2.5° around declination =0 can 
be found at http : / /www. sdss3.org/gallery/sdssjpie2.php . Other surveys 
have the following opening angle in declination : 3° for the second Center 
for Astrophysics (CFA2) survey, see Geller & Huchra (1989) , 4° for the 



2dF Galaxy Redshift Survey (2dFGRS), see Colless et al 



the Six-degree Field instrument (6dF), see Jones et al. (2004). Particular 



(2001) and 4° for 



attention should be paid to the fact that the astronomical slices are not a 
plane which intersects a Voronoi Network. 

In order to quantify this effect we introduce a confusion distance, D c , as 
the distance after which the half altitude of the slices equalizes the observed 
average diameter D obs 

(9) 



D, 



tau(-) 



_JJobs 

2 



where p is the opening angle of the slice and D obs the averaged diameter of 
the voids. In the case of SDSS tp = 2.5° and therefore D c = 835 .4/ h Mpc 
when D obs = Mpc. For distances greater than D c the voids in the 
distribution of galaxies are dominated by the confusion. For distances lower 
than D c the filaments of galaxies can be considered the intersection between 
a plane and the faces of the Voronoi Polyhedrons. The PVT offers a new 
classification for extragalactic aggregates. 

1. Groups and clusters of galaxies are situated on the faces of the PVT 
with typical linear dimensions of ~ 16 Mpc. 

2. The local super-cluster is made by the network of many faces of the 
PVT. For distances equal to or slightly greater than D obs /2, a spherical 
cut should reveal 16 voids . 

3. The filaments of galaxies visible in the slices of oriented catalogs are 
due to the intersection between a plane and the PVT network of faces 
as first approximation. An improvement can be obtained by coding the 
intersection between the slice of a given opening angle and the PVT 
network of faces 



sec 



Zaninetti (2006, 2010). As an example Figure 3 



reports both the CFA2 slice as well the simulated slice. 



4. The Great Wall, see Geller & Huchra (1989), is due to a combined 



effect of the photometric maximum, see Equation (19) and the PVT 
network of faces. 
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Figure 3: Polar plot of the real galaxies (green points) belonging to the sec- 
ond CFA2 redshift catalog and the simulated galaxies in the PVT framework 
(red points). More details can be found in Zaninetti (2006). 



The cellular structure in the spatial distribution of galaxies can be gen- 
erated in other approaches such as the long-range force acting only between 



non-baryonic particles, see Keselman et al. (2010). 



3.1 Thickness 

The concept of thickness of a face can be introduced by a numerical experi- 
ment. We generate N r random galaxies in 2D and then we progressively shift 
toward 3D increasing the thickness of the box from zero to a given value. The 
analysis of the number of galaxies, N(R), as function of the distance can be 
done by assuming a power law dependence of the type 



N(R) = CR? 



(10) 



The two parameters C and /3 can be found from the following logarithmic 
transformation 

\n(N(R)) = ln(C) + (3 \n(R) , (11) 
which can be written as 

y = a LS + b LS x . (12) 
The application of the least squares method through the FORTRAN subrou- 



tine LFIT from Press et al. (1992) allows of finding a^s, b^s an d the errors 



a a and o^. The resulting dimension increases progressively with thickness, 
see Figure |4| 
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thickness (Mpc) 



Figure 4: The dimension of the number of galaxies represented with vertical 
error bars as function of the thickness in Mpc. 



The fact that the dimension is not exactly 2 at thickness =0 , more exactly 
we have /3 = 1.8 when A/" r =10000 , is due to the cardinality of N r , as an 
example we have /3 = 1.96 when iV r =1000000. This first analysis indicates 
that the dimension is function both of a geometrical factor as thickness and 
numerical factor , i.e. the number of considered galaxies. We now outline the 



method that allow us to insert the galaxies on the thick faces, see Zaninetti 



(1995, 2010). A practical implementation is to model a decreasing probability 



of having a galaxy in the direction perpendicular to the face. As an example 
we assume a probability, p(x), of having a galaxy outside the face distributed 
as a Normal (Gaussian) distribution 



p(x) 



X 



a ( 27r )i/2 ex P- 2a 2 

where x is the distance in Mpc from the face and a the standard deviation in 
Mpc. Once the complex 3D behavior of the faces of the Voronoi Polyhedron 
is set up we can memorize such a probability on a 3D grid P(i,j,k) which 
can be found in the following way 

• In each lattice point (i, j, k) we search for the nearest element belonging 
to a Voronoi face. The probability of having a galaxy is therefore 



(13) 



computed according to formula (13). 



A number of galaxies, — n* x side 3 is then inserted in the box; here 
n* represents the density of galaxies 
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Figure 5: 3D visualization of the edges of the Poissonian Voronoi-diagram. 
The parameters are pixels= 60, N s = 60, side = 89 Mpc, h = 0.7 and 
amplify= 1.2. 

Figure [5] visualizes the edges belonging to the Voronoi diagrams and Figure [6] 
represents a cut in the middle of the probability, P(i,j, k), of having a galaxy 
to a given distance from a face. 

We are now ready to evaluate the spatial dimension in two cases 

• ntrials centers randomly situated on a face , 

• ntrials centers randomly situated on a edge, 

see Table [TJ The case of evaluation starting from the edges produces a bigger 
dimension because we enclose three faces rather than one. 

4 Observations 

This section processes the Sloan Digital Sky Survey Data Release 7 (SDSS 
DR7) and the Two-Micron All Sky Survey (2MASS). 
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minim um = 



0,0000 



m o x i m u m = 



0.1995 




Figure 6: Cut in the middle of the 3D grid P(i,j,k) which represents 
a theoretical 2D map of the probability of having a galaxy. The Voronoi 
parameters are the same as in Figure [5] and a = 0.8Mpc. The X and Y units 
are in Mpc. 



Table 1: The dimension in two cases when ntrial=10. 



case 


value 


faces 


1.65 ±0.1 


edges 


1.79 ±0.07 
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Table 2: The statistical parameters of the effective radius in SDSS DR7. 



parameter value 

mean 18.23/i -1 Mpc 

variance 23.32/i -2 Mpc 2 

standard deviation 4.82/iT 1 Mpc 

kurtosis 0.038 

skewness 0.51 

maximum value 34.12/i _1 Mpc 

minimum value 9.9h~ l Mpc 



4.1 Observed Statistics of the Voids 

The distribution of the effective radius between the galaxies of SDSS DR7 



has been reported in Pan et al. (2012) and is also possible to extract the 



catalog which is visible at 

http : //www. physics . drexel . edu/~pan/voidcatalog.html 

; this catalog contains 1054 voids and Table [2] reports their basic statistical 
parameters. 



4.2 The 2MASS 

The 2MASS is a catalog of galaxies which has instruments in the near-infrared 
J, H and K-bands (1-2.2 \i m) and therefore detects the galaxies in the so 
called "Zone of Avoidance," see iJarrettl (120041); ICrook et all (120071). At the 



moment of writing the 2MASS Redshift Survey (2MRS) consists of 44599 



galaxies with redshift in the interval < z < 0.09, see Huchra et al. (2012). 



Figure [7] reports a spherical cut at a given radius in z of the LSC accord- 
ing to 2MRS, which is available online at 

https://www.cfa.harvard.edu/ huchra/seminar/lsc/lsc.dat . The LSC is cen- 
tred on the Virgo cluster , has a flattened structure and it's radius is ~ 40 
Mpc. 
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Zrf„=0.a 133+1 3^= 0.014404 No galaxies = 693 




Figure 7: Hammer-Aitoff projection in galactic coordinates of a spherical 
cut of the LSC data at 0.0133 < z < 0.0144 or 5G.9Mpc <D< 61.67 . 

A way to smooth the results of Figure [7] is through contour of isodensity 
of galaxies in a cube of sides equal to 16 Mpc at a given distance, see Figure |8j 
The previous figure contains voids in which the density of galaxies is zero. 
In order to have in all regions of the spherical cut a minimum density of 1 
galaxy/ (chosen volume) we consider the density as evaluated in a cube which 
has a side equal the averaged diameter as given by SDSS DR7, see Figure [9j 
Now the density of galaxies has a minimum value of 1 at the center of a void 
and a maximum value of 270 on a side which comes out from the intersection 
of a sphere with a PVT. This means that the density of galaxies varies over 
more than two decades on a spherical cut. 

The concept of homogeneity can be associated to a random sample which 
has the same cardinality of the astronomical galaxies but placed at random 



in the same analyzed shell; Figure 10 reports an example of such random 
galaxies. The ratio 

H(r) = No{r) (14) 

quantifies the homogeneity of the astronomical sample with Nc{r) represent- 
ing the number of real galaxies and A#(r) the random galaxies in a cube of 
side 2r , (Nf>(r)) represents the number of random galaxies averaged over 
91 x 181 directions. A value of of if = 1 means that the sample is homo- 
geneous , a value lower that 1 means that the sample is not homogeneous. 
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Figure 8: Logarithmic contours of density of galaxies of the Large Super 
Cluster. The density is evaluated in cubes of side = 16 Mpc at z — 0.0133 
or D = 56.9Mpc. 



Figure 9: Linear contours of density of galaxies of the Large Super Cluster. 
The density is evaluated in cubes of side = 16 Mpc at z — 0.0133 or D = 
56.9Mpc. 
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Figure 10: Hammer-Aitoff projection in galactic coordinates of a spherical 
cut of a random sample of galaxies at 0.0133 < z < 0.0144 or 56.9Mpc < 
D < 61.67 . The number of galaxies is the same of Figure [7} 



An example of the various value of H is reported in Figure 10 and Figure 



12 reports the numerical survival function for the grid points. From this last 
figure is possible to conclude that we have non homogeneity over the 68 % 
of the grid points of the spherical cut. 

Figure [9] which shows the linear density contours of galaxies and Figure 
[TT1 which shows the linear contours of H are similar because the numbers of 
simulated galaxies in the selected cube of the random sample fluctuates less 
in respect to the sample of galaxies. 

5 Luminosity effects 

This section reports the luminosity function for galaxies, the Malmquist bias 
which allows to build a complete sample of galaxies, and the numerical de- 
termination of the spatial dimension of the LSC 

5.1 Luminosity Function 

A model for the luminosity of galaxies is the Schechter function 

mdL=(~)(~TeM-~)dL , (15) 
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galactic 



Figure 11: Linear contours of H as given by (14) in the Large Super Cluster. 
The density of the astronomical and random samples are evaluated in cubes 
of side = 16 Mpc at z = 0.0133 or D = 5Q.9Mpc. 




Figure 12: Survival function of of the ensemble of 91 x 181 grid points as 
function of H at D = 5Q.9Mpc. The median is at H— 0.6 , in an homogeneous 
sample it should be at H— 1. 
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Table 3: The parameters of the Schechter function and bolometric magni- 
tude for the 2MRS in the K s — band. 



parameter 2MRS 

M* - 5 log 10 h [mags] (-23.44 ± 0.03) 

a -0.96 ±0.05 

$* [h 3 Mpc- 3 ] ((1.08±0.16)10~ 2 ) 

M% s 3.39 

h 0.7 



where a sets the slope for low values of L, L* is the characteristic luminosity, 
and $* is a normalization. This function was suggested by Schechter (1976) 
and the distribution in absolute magnitude is 



<&(M)dM 



(0.4Znl0)$*10 a4(a+1)(M *- M) 
x exp(-10 a4(M *- M) )rfM , 



(16) 



where M* is the characteristic magnitude as derived from the data. The 
parameters of the Schechter function concerning the 2MRS can be found in 



Cole et al. (2001 ) and are reported in Table 3 



The number of galaxies at a given flux / as a function of the redshift z, 



see Formula (1.104) in Padmanabhan (1996), is 



dN 
dQdzdf 



4tt( 



Hn 



^crit 



(17) 



where dQ, dz, and df are the differentials of the solid angle, the red-shift, 
and the flux, respectively, and 



^crit 



HqL* 
Airfc 2 



(18) 



where c is the velocity of light. The number of galaxies at a given flux has a 



maximum at z 



where 



ZcritVa + 2 



(19) 
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Figure 13: The galaxies of the 2MRS with 7.23 < m < 8.92 or 
4881116 M ^ 2 < / < 23044324 M ^ 2 are organized in frequencies versus helio- 
centric redshift, (empty circles); the error bar is given by the square root of 
the frequency. The maximum frequency of observed galaxies is at z = 0.0058. 



The full line is the theoretical curve generated by (17), the Schechter function 
of luminosity, with parameters as in Table [3] In this plot Mq = 3.39 ,h = 
0.7 , x 2 = 132.93 and the number of bins 15. 



Figure [13] reports the number of observed galaxies in the 2MRS catalog at a 
given apparent magnitude and the theoretical curve as represented by Equa- 
tion (17). The merit function \ 2 can be computed as 



^2 _ •y^^ n theo( z ) n astr{ z ) ^Q) 
j=l ^riastriz) 

where n is number of data , the two indexes theo and astr stand for the- 
oretical and astronomical, respectively and <J nast r(z) 2 * s the variance of the 
astronomical number of data; the obtained value is reported in the caption 



of Figure 13 



On analyzing Figure 13 we suggest that the disagreement between theory 
and observations is due to the fact that the number of galaxies does not in- 
creases exactly with R 3 as inserted in the theory but with R 1S7 , see equation 



(26), as given by the cellular structure. 
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Figure 14: The absolute magnitude M of 5512 galaxies belonging to the 
2MRS when M & = 3.39 and H = 70 km s" 1 Mpc -1 (green points). The 



upper theoretical curve as represented by Equation (21) is reported as the 
red thick line when m£=11.19. 



5.2 Malmquist Bias 



The Malmquist bias, see Malmquist (1920, 1922 



to the stars and later on to the galaxies by Behr 



, was originally applied 
(1951). The observable 



absolute magnitude as a function of the limiting apparent magnitude, m L , is 



M L = m L - 5 log u 



c z 
Ho 



25 



(21) 



The previous formula predicts, from a theoretical point of view, an upper 
limit on the maximum absolute magnitude which can be observed in a catalog 
of galaxies characterized by a given limiting magnitude and Figure |T4"| reports 
such a curve and the galaxies of the 2MRS. 

The interval covered by the LF of galaxies, AM, is defined by 



AM = M max - M mm 



(22) 



where M max and M min are the maximum and minimum absolute magnitude 
of the LF for the considered catalog. The real observable interval in absolute 
magnitude, AMl, is 



AM = Mr 



Mr, 



(23) 
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We can therefore introduce the range of observable absolute maximum mag- 
nitude expressed in percent, e a (z), as 

€s{z) = w x 100% • (24) 

This is a number which represents the completeness of the sample and, given 
the fact that the limiting magnitude of the 2MRS is m£=11.19, it is possible 
to conclude that the 2MRS is complete for z < 0.00016 . This efficiency 
expressed as a percentage can be considered a version of the Malmquist bias. 

The previous treatment does not take into account the fT-corrections 
between observed (Q) and desired (R) bandpasses, Kq R (z) , see equation 



(1) Blanton & Roweis (2007). The i^-correction depends from the chosen 
catalog and bandpass . In our case we are interested in band K s in the 
2MASS and the following approximate relationship can be drawn from Figure 



14 in Blanton & Roweis (2007): 



K K {z) pa 1.22 . (25) 

As an example at z ~ 0.016 , which is the maximum processed value of 
redshift of the LSC, the /^-corrections are 0.02 mag. This means that the 
.ff -corrections are negligible when applied to the LSC , which is a subset of 
the 2MASS. 



5.3 The Radial Distribution of Galaxies 

In a 3D random distribution, the number of galaxies grows oc R 3 where R is 
the distance from the center (our galaxy) in Mpc. In a framework of galaxies 
situated on the first face of a PVT their number should grow oc R 2 assuming 
that R denotes the distance from the center of the face. The counting of 
the nearest galaxies can be done in light of the Malmquist bias and the first 
face of a PVT. According to (J8J), the galaxies on the first face will terminate 
at Dp pa 8.75 Mpc which corresponds to z — 0.00204. The Malmquist line 



as represented by (21) suggests that the sample of galaxies in the interval 



< z < 0.00204 and characterized M < —18.51 is complete. According to 



equation (10) the spatial behavior is 

N(R) = (3.52 ± OA9)R^ 87±om ^ (26) 
when < z < 0.00204 and M < -18.51 , 
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Figure 15: The number of galaxies as function of the distance in Mpc for the 
2MRS catalog (full line) and observed number of galaxies with vertical error 
bars. The error is computed as the square root of the number of galaxies in 
the considered bin. 



and Figure 15 reports such a relationship. 

This last result can be tested on a random 2D sample in which the galaxies 
are generated on a plane having the same area of the astronomical analyzed 
region. The result of the simulation for the random galaxies gives 

N(R) = (1.43 ± O.OS)^ 2 - 03 ^- 082 ) (27) 
when < z < 0.00204 . 



A comparison of the two curves of growth ( random galaxies and observed 
galaxies) is reported in Figure 16 

Others approaches use the generalized fractal dimension estimator to (i) 
analyze the scaling properties of the clusters of galaxies, see Borgani et al.| 
(1993), (ii) determine the scale of homogeneity in the Universe as derived 
from various galaxy redshift surveys, see Bagla et al. (2008), (iii) study the 
transition to homogeneity, defined as the scale above which the fractal dimen- 
sion of the underlying galaxies distribution is equal to the ambient dimension 
of the space in which the galaxies are distributed, see Yadav et al. (2010). 
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Figure 16: The number of real galaxies as function of the distance in Mpc 
for the 2MRS catalog (full line) and number of random galaxies generated 
on a plane (dashed line). 

6 Conclusions 



In Section we learned that the concept of homogeneity does not corre- 
spond to the observations of LSC because the 63 % of 91 x 181 grid points 
processed at a distance of w 59 Mpc have the indicator of homogeneity lower 
than 1. The concept of a non-isotropic universe is a consequence of the pre- 
vious statement, in other words, the appearance of the universe depends on 
the chosen observer /galaxy. These facts, based on the theoretical and ob- 
servational bases, allow us to assert a local new version of the Cosmological 
Principle "The galaxies, as well our galaxy, are situated on the faces of a 
Voronoi polyhedron. The spatial distribution of galaxies will follow approxi- 
mately the geometrical rules as well the standard photometric rules." A first 
set of consequences can be drawn 

1. The number of galaxies in the first face of a PVT should grown as 
N(R) oc R 2 . An analysis of the number of galaxies in the first 8.75 
Mpc of the LSC gives N(R) oc R 187 . 

2. A new definition of the density of galaxies on the faces of a PVT, rif, 
which in the framework of thick faces is 

nf « n— , (28) 
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where p = s/R is a dimensional parameter lower than 1, s is the thick- 
ness of the faces, and n is the usual density of galaxies. 



3. The so called "filaments of galaxies" are here classified as 2D sections 
of a 3D PVT which are made by irregular polygons. 

4. The reformulation here presented is in agreement with the anisotropies 
in the short and intermediate gamma ray bursts which are explained 
by a breakdown of the Cosmological Principle, see Meszaros et al. 



(2009),Meszaros et al. (2009). 



5. The theoretical dimension which regulates the number of galaxies as 
function of the distance is different when the center of the selected 
sphere belongs to a thick face or a thick edge. 

Different approaches which analyze the SDSS seem to suggest the validity 
of cosmological principle and they are: (i) the test the homogeneity of the 
universe at z « 0.3 with the luminous red galaxy (LRG) spectroscopic sample 
of the SDSS, see Hogg et al. (2005), (ii) a multifractal analysis of the galaxy 



distribution in a volume-limited sub-sample from the (SDSS) which finds that 
the galaxy distribution becomes homogeneous at a length-scale between 60/h 



Mpc and 70/h Mpc, see Sarkar et al. (2009). We remember that 70/h Mpc as 



a linear dimension comprises ~ 2 PVT cells. The analysis here presented is 
limited to galaxies within 8 Mpc in order to measure the spatial dimension. 
Tests carried out by other authors also do not find homogeneity within 8 
Mpc. Rather, the transition to homogeneity is found to occur at a much 
larger length-scale. So there does not appear to be any inconsistency with 
the findings of our paper and the currently accepted picture of the transition 
to homogeneity. In the framework of the cellular universe the transition from 
non homogeneous to homogeneous regime is reached when the cube in which 
the density of galaxies is evaluated comprises many cells. 
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